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Unambiguous discrimination of mixed quantum states: 
optimal solution and case study 
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We present a generic study of unambiguous discrimination between two mixed quantum states. We 
derive operational optimality conditions and show that the optimal measurements can be classified 
according to their rank. In Hilbert space dimensions smaller or equal to five this leads to the complete 
optimal solution. We demonstrate our method with a physical example, namely the unambiguous 
comparison of n quantum states, and find the optimal success probability. 
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According to the laws of quantum mechanics, two non- 
orthogonal quantum states cannot be distinguished per- 
fectly. This fact has far-reaching consequences in quan- 
tum information processing, e.g. it allows to generate a 
secret random key in quantum cryptography. In spite of 
the fundamental nature of the problem of state discrimi- 
nation, determining the optimal measurement to distin- 
guish two (mixed) quantum states is far from being triv- 
ial. 

In the literature, two main paths to state discrimina- 
tion have been taken [l|: firstly, in minimum error dis- 
crimination, the unavoidable error in distinguishing two 
states from each other is minimized. This problem has 
been completely solved in Ref. |2J. Secondly, in unam- 
biguous state discrimination (USD), no error is allowed, 
but an inconclusive answer may occur. The optimal USD 
measurement minimizes the probability of an inconclu- 
sive answer [J-Q. Although USD has found much atten- 
tion in the recent years, and special examples have been 
solved, no general solution is known so far for the case of 
mixed states. A strategy that is analogous to USD, but 
applicable also for linearly dependent states are maxi- 
mum confidence measurements, discussed in Ref. [6[. 

The aim of the current contribution is to present the 
optimal USD measurement for cases which cannot be re- 
duced to the pure state case and thus to known solu- 
tions. This analysis can be applied for the unambiguous 
discrimination of any two density operators acting on a 
Hilbert space up to dimensions five. This goes beyond 
previous results which require a high symmetry or other 
very special properties of the given states |7Hl3l] . We will 
show the main ideas and steps towards the solution, and 
explain the technical details elsewhere [14| . 

The scenario of optimal unambiguous discrimination of 
two density operators is as follows: two (normalized) den- 
sity operators g\ and £>2, acting on a finite-dimensional 
Hilbert space J^f occur with a priori probability p\ and 
P2, respectively, where pi + P2 — 1- We will denote 
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the support of a density operator g as the orthocom- 
plement of its kernel, (supp g) = ker g. A measure- 
ment for USD is described by a positive operator valued 
measure (POVM), i.e., a family of positive semi-definite 
operators {E\,E2,E->} with E\ + E% + E> = 1, obey- 
ing the constraints for unambiguity, Xx(EiQ\) = and 
Xx{E\Q>i) — 0. The operator E? corresponds to the in- 
conclusive outcome while E\ and E^ correspond to the 
successful detection of g\ and 02, respectively. The aim is 
to find a POVM which maximizes the success probability 
-Psucc = Pi tr(EiQi) + P2 tr(£ , 2g2)- Let us introduce here 
the useful notation 71 = pigi and 72 = P2Q2- Thus, the 
success probability reads P SUC c = tr(Ei'ji) + ti(E 2 ^2)- 

What are the relevant structures of the density op- 
erators and measurement operators? The unambiguity 
condition tr(i?27i) = means that the support of E2 is 
a subspace of the kernel of 71 . The second unambiguity 
condition reads supp£a C ker 72. Obeying these con- 
straints, one has to maximize the sum of the scalar prod- 
ucts tr(i?i7i) and tr (^272)5 while keeping E? positive. 
Due to the reduction theorems in Ref. [8(, the optimiza- 
tion problem reduces to the case of a strictly skew pair of 
(unnormalizcd) density operators. The operators 71 and 
72 are called strictly skew, when they neither possess any 
parallel component, i.e., supp 71 fl supp 72 = {0}, nor any 
orthogonal components, i.e., supp7i nker7 2 = {0} and 
supp 72 fl ker7i = {0}. A simple example for a strictly 
skew pair of unnormalized density operators is any pair 
of pure states, 71 = p|0iX0i| and 72 = (1 - p)|02)(02|, 
with < I (0i |02) I < 1 and < p < 1. Both opera- 
tors of such a strictly skew pair have the same rank, and 
the sum of both ranks cannot exceed the dimension of 
the underlying Hilbert space. - - Below we will show a 
constructive method to discriminate two skew density op- 
erators of rank two. This solves optimal USD in all cases 
where one of the given states has rank two, and hence 
in articular the case with a Hilbert space of dimension 
smaller equal five. 

In the following we will only consider skew pairs of 
unnormalized density operators and proper USD mea- 
surements. We call a USD measurement proper, if it sat- 
isfies supp(i?i + E 2 ) C supp (71 +72). It is sufficient to 



only consider proper measurements, since the subspace 
ker 71 n ker 72 cannot contribute to the success probabil- 
ity Q. 

In Rcf. [16] Eldar and collaborators showed that the 
optimality of a USD measurement can be proved via the 
existence of a certain operator that fulfills a set of condi- 
tions. However, no constructive way to find this operator 
was provided. Starting from these conditions we derive 
the following set of necessary and sufficient requirements 
for the optimality of a proper USD measurement: 



£?(72~7i)£?(l-£?) = 0, 
Ai£ ? (72-7i)£?A 2 = 0, 
Ai-E?(7 2 - 7i)-E?Ai > 0, 
A 2 £?(7i -li)Ei A 2 > 0. 



(la) 
(lb) 
(lc) 
(Id) 



Here, Ai is the projector onto ker 72, and A2 is the pro- 
jector onto ker 7!. The details of the derivation are pre- 
sented elsewhere [14( . Note that the methods used in or- 
der to arrive at Eqns. (JTJ cannot be generalized to the dis- 
crimination of more than two states. (For special cases, 
however, cf. Ref. (I7J.) 

Let us point out two observations from Eqns. ([1]): first, 
neither E\ nor E%, but only the operator E? appears in 
this set of equations. This is due to the fact, that from 
E? it is possible to uniquely reconstruct E\ and E 2 , as 
Ei'ji — 7,: — E-r'fi holds for i = 1,2. Second: neither 
7i — 72 > nor 72 — 71 > can hold for a strictly 
skew pair of operators, and thus it is non-trivial to fulfill 
Eq. (JTcJ and Eq. ([Td]). The set of equations (fTa| - (fTd"|) 
provides an efficient tool in optimal USD: one might 
be able to guess a measurement, e.g. from the symme- 
try of a given USD problem, and can then verify easily 
whether it is optimal. Moreover, one can use these equa- 
tions in a constructive way in order to find the solution 
for Ei, which then uniquely defines an optimal POVM. 
Below, we will show explicitly how to construct the opti- 
mal measurement from Eqns. ([1]) for the example of state 
comparison. 

It has been an open question whether the optimal 
USD measurement is unique. This is indeed the case. 
The structure of the proof is as follows: As pointed out 
above, a USD measurement is already defined via E?. It 
can be shown |14| that for optimal proper measurements 
the rank of E? is fixed, namely rankE? = rank(7i72) + 
dim ker (71 +72). Assuming that there would be two 
optimal operators E? and E,, their convex combination 
|(-E? +E' ? ) would also describe an optimal measurement. 
However, for positive semi-definite operators E? and E!,, 
the identity rank(E? + Elf) = rankE? = rankE? can 
only hold if suppE? = suppE-'. When the support of 
E? is given, the operator E? is uniquely determined via 
Eq. ([lap . Thus, the optimal proper USD measurement is 
unique. 

The uniqueness of the optimal measurement now al- 
lows a meaningful characterization of the optimal USD 
measurement. We introduce a classification of the differ- 
ent types of optimal USD measurements according to the 



rank of the measurement operators E\ and E2 ■ A mea- 
surement type is specified by (rank E±, rank E2). This 
classification turns out to be vital for the construction 
of optimal measurement strategies from Eqns. ([T}. For 
given density operators Q\ and Q2 and a given a priory 
probability p\ — 1— P2, one particular measurement type 
is optimal, due to the uniqueness of the optimal solu- 
tion. While varying p\ some or all of these measurement 
types may occur, see Fig. [1] for an illustration. With 
r = rank 71 = rank 72, one arrives at the constraints 
rank .Ei < r, rankE2 < r, and 



r < rank E\ + rank E 2 < 2r. 



(2) 



Eq. @ follows from the geometry of unambiguous mea- 
surements and the fact that in the optimal case rank E? = 
dim ker 7172 holds. The two extremal cases where either 
the lower or the upper bound in Eq. (J2J) is reached cor- 
respond to special situations. 

The case of the upper bound in Eq. ^j, where 
rank Ei = r = rankE2, is the well- understood fidelity 
form measurement: Intuition might tell that the suc- 
cess probability should be a function of some distance 
measure between the two states (this is indeed true 
for minimum error discrimination, where the smallest 
achievable error probability is a function of the trace 
distance between the unnormalizcd density operators). 
Here, for the case with rankEi = r = rankE2 the suc- 
cess probability is the square of the Bures distance, i.e., 
Pfid = l-2tr|V7iV72| IS EL HE! (while, in general, 
Pfid is an upper bound on the success probability |15|). 
In fact, formally the construction of the fidelity form 
measurement is always possible [1JJ and the resulting 
operator E? always satisfies all conditions in Eqns. ((T|). 
However, this operator in general fails to satisfy the con- 
dition 1 — E? > 0. The measurement types for which 
rank E\ + rank E2 < 2r occur due to this very positivity 
condition. In a geometric language the optimal measure- 
ment is on the border of the allowed (positive) measure- 
ments, unless rank Ei = r = rankE2. One can compute 
two numbers p\ ow and p up for given q\ and Q2 , such that 
the fidelity form measurement is optimal if and only if 

Plow <Pl <Pup- 

In the case of the lower bound of Eq. ©, where 
rank Ei + rankE2 = r, the operators E\, E2, and E? 
are projectors, i.e., the optimal measurement is a von- 
Neumann measurement. A special situation occurs when 
rank Ei = and rankE2 = r or rankEi = r and 
rankE2 = 0. This is interpreted as follows: For very 
small p\ it will turn out to be advantageous to ignore 
Q\ by choosing E\ = 0. This case is referred to as sin- 
gle state detection of Q2, because the state Q\ is never 
detected. As then E? = 1 — E 2 , from Eqns. ([T]) only 
Eq. (|lc|) remains, and this inequality can be written as 



71(72 -71)71 > 0. 



(3) 



The success probability for single state detection of 72 is 
given by P succ = tr(A272), where A2 was defined above 




FIG. 1: USD measurement types for r — 2, as allowed by 
the constraint in Eq. (|5J. Projective measurements are indi- 
cated by squares, non-projective ones by circles. The arrows 
illustrate an example for a possible path between the measure- 
ment types, while the probability p2 is varied from P2 = to 
P2 = 1. The start point is necessarily type (2, 0) and the end 
point type is (0, 2). The types (0, 2), (2, 0), and (2, 2) will only 
occur once. Which other types are visited inbetween, and in 
which order, depends on the concrete example. 



as the projector onto ker7i. Eq. ([3]) implicitly defines 
a calculable threshold for pi, below which it is advanta- 
geous not to detect g±. This threshold is always larger 
than 0, i.e., single state detection is always optimal for 
a finite regime. Analogous considerations hold for small 

P2- 

So far our considerations have been independent of r. 
Let us now consider specific values for r. For r = 1, 
i.e., the case of pure states, only the single state detec- 
tion measurement or the fidelity form measurement may 
occur. Hence the problem of unambiguous discrimina- 
tion of pure states is well understood 18]. Furthermore, 
any USD task where the two density operators can si- 
multaneously be brought in a diagonal form with 2x2- 
dimensional blocks (the "block-diagonal" case), can also 
be solved by treating the corresponding orthogonal sub- 
spaces independently [3, [HI, LL9| ■ For all other cases only 
solutions for special cases are known [irj| - [l3{ . - For 
r = 2 there are six possible measurement types which 
are summarized in Table |U The optimal measurements 
for the types (1,2), (2, 1), and (1,1) remain to be deter- 
mined. For each of these types, the Eqns. ([T]) reduce to a 
polynomial equation [lj] and hence the analytic solution 
for the case r = 2 is completed. 

Let us now study the important example of quantum 
state comparison and demonstrate explicitly how to solve 
Eqns. (Q} for the case of measurement type (1,1) which 
occurs for a wide range of parameters. We consider state 
comparison of n pure quantum states, where each of the 
states is taken from the set {|'0i)> IV^)}, with correspond- 
ing a priori probabilities {771 , 772}, m + 772 = 1. In quan- 
tum state comparison [lCl Il5l I19M22J one aims at answer- 
ing the question whether the given n quantum states are 



rank E\ 


rank E2 


type 





2 


(0,2) 


1 


2 


(1,2) 


2 


2 


(2,2) 


1 


1 


(1,1) 


2 


1 


(2,1) 


2 





(2,0) 



properties 



single state detection, projective 
non-projective measurement 
fidelity form measurement, non-pro j. 
projective measurement, cf. example 
non-projective measurement 
single state detection, projective 



TABLE I: Measurement types for the case r 
about the properties see main text. 



2. For details 



equal or not. Applications of this task in quantum infor- 
mation are e.g. quantum fingerprinting 23| and quantum 
digital signatures [2J] . For n = 2 the optimal unambigu- 
ous measurement for quantum state comparison has been 
given in Ref. [10|, [22|- For n > 3, the corresponding USD 
task reduces to the unambiguous discrimination of two 
mixed states of rank 2, i.e., r = 2. 

State comparison of n states is equivalent to the dis- 
crimination of (cf. Ref. [221 ]]) 



\0n 



i c = {vMi)^i\r n + {m\^)m) 

7d = (r?i|MVi|+%hW2|)®"-7. 



(4) 
(5) 



Due to Theorem 2 in Ref. [8fl it remains to consider the 
reduced operators 7^ and 7^ , which are given by the pro- 
jection of 7e and 7d onto (supp 7 C + ker7d), respectively. 
It is straightforward to see that for n > 3 this discrimi- 
nation task cannot be reduced further and that no block- 
diagonal structure is present unless 771 = 772 = \. 



We next construct a basis of supp 7 e and of ker 7d . 
convenient basis of supp 7 C is given by 



A 



\Qjn 



\4> ± ) ex \^r n ± i^> 

We define c = (ipilipz) with < c < 1. Using \ipi 
\ip2) ~ <#i) and \ip£) ex \ipi) - c\ip 2 
be constructed as 



(6) 



a basis of ker 7d can 



IwiJoclViT ±l^> 



<gm 



(7) 



Gram-Schmidt 
w_)} yields 



orthogonalization of 
the orthonormal ba- 
+ ker7 d . 
basis of 
is an or- 
In fact, 
of these 



Now c 

{|0+),|<M> + 

sis {|0 + ),|0_),|er + ),|cr_)} of supp7 
Then {|cr + ), |cr_)} is an orthonormal 
ker 7 ^ n supp(7 c r + 75) while {|w+},|cj_}} 
thonormal basis of ker7^ n supp(7^ + 7^} 
they form Jordan bases (cf. e.g. Ref. [15|, [2j_.,, 
subspaces, i.e., (<j t \uj±) — 0. The remaining overlaps 
(<r±|u;±) are equal for odd n [degenerate Jordan angles). 
We now study for general but odd n > 3 the solution 
of the conditions in Eqns. ([T]) while restricting our 
considerations to the measurement type (1, 1). 

The measurements of type (1,1) are von-Ncumann 
measurements, where E e and E& both have rank 1, i.e., 
E c = |XeXXe| and E d = |xd)<Xd|- In particular the vec- 
tors |xe) and |xd) must be orthogonal and normalized. 
We use the parametrization \x c ) oc |w+) + x*\lo^) and 
\\d) oc x\a + ) — |<7_) where a; is a complex variable [271 ] . 




FIG. 2: Maximal success probability for comparison of 3 pure 
quantum states, taken from the set {|^>i), l^)}, as a function 
of the a priori probability ?7i and the overlap {ipi\ij)2). Darker 
areas correspond to lower success probability. The dashed line 
indicates the bound from the conditions (|lc|) and (|ld[) . 



We now evaluate the necessary and sufficient condi- 
tions for optimality in Eqns. (p}: Eq. (fla)) is satisfied 
for any x. Let us abbreviate (w |7 e |w(,) = G° b and 
(cra|7d|crfc) = Gf, where a, b E {+,-}. Eq. ((TEJ) now 
becomes a scalar equation which is only quadratic in x; 
in matrix notation Eq. (jlb[) reads 



(l,x)(G c -G d )(-a ; ,l) T = 0. 



(8) 



Similarly, the positivity conditions (|lc|) and (jld[) simplify 
to scalar inequalities. 

With the help of a computer algebra system we obtain 
for n = 3 the optimal success probability 



the dashed line the optimal measurement is of type (1,1) 
and the success probability is given by Eq. ([9]). We find 
from numerical analysis that the optimal measurement 
is a fidelity form measurement in the remaining cases. 
Note, that for a wide range of the parameters the optimal 
measurement is a von-Neumann measurement and hence 
may be implemented physically without the need of an 
auxiliary system. 

In summary, we have presented a strategy to find the 
optimal measurement for unambiguous discrimination of 
two mixed quantum states acting on a five-dimensional 
Hilbert space. Our method can in principle also be ap- 
plied to the discrimination of two quantum states in gen- 
eral dimensions. Our results are also useful in other con- 
texts, e.g. quantum state filtering: in Ref. Qj it has been 
shown how to optimally distinguish between one pure 
state from a given set and the remaining ones. With our 
method one could filter a subset of states from the whole 
set. In connection to quantum algorithms, one could thus 
distinguish between two sets of Boolean functions, rather 
than between one function and a set of functions. The 
results presented in this paper could also be used to prove 
optimality for the universal programmable state discrim- 
inator suggested in Ref. [2g. As the optimal measure- 
ment is unique, the optimal device discussed in Ref. [26[ 
cannot be simplified. Furthermore, in Ref. |11| the im- 
portance of unambiguous discrimination in the context of 
quantum key distribution was shown with particular em- 
phasis on the case of states of rank two. As an outlook, 
our strategy seems a promising path for the generaliza- 
tion to unambiguous state discrimination of more than 
two states. 
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F4c 2 + l)a+ 
+ (l-c 2 )(2 + 



W)}, (9) 



withW = [(1-c 6 



-4(1 



')](!- a) + a 2 c 2 



a = 4 771 772 ■ Note that this expression is only valid if in 
addition the inequalities (fTc)) and (lid)) hold. The success 
probability is illustrated as contour plot in Fig. [2] Above 
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